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Abstract 

In this paper, we report the results of constraining the holographic dark energy model with spatial curvature and 
massive neutrinos, based on a Markov Chain Monte Carlo global fit technique. The cosmic observational data include 
the full WMAP 7-yr temperature and polarization data, the type la supernova data from Union2.1 sample, the baryon 
acoustic oscillation data from SDSS DR7 and WiggleZ Dark Energy Survey, and the latest measurements of Hq from 
HST. To deal with the perturbations of dark energy, we adopt the parameterized post-Friedmann method. We find that, 
for the simplest holographic dark energy model without spatial curvature and massive neutrinos, the phenomenological 
parameter c < 1 at more than Act confidence level. The inclusion of spatial curvature enlarges the error bars and leads 
to c < 1 only in about 2.5<x range; in contrast, the inclusion of massive neutrinos does not have significant influence 
on c. We also find that, for the holographic dark energy model with spatial curvature but without massive neutrinos, 
the 3cr error bars of the current fractional curvature density Q^o are still in order of 10" 2 ; for the model with massive 
neutrinos but without spatial curvature, the 2<x upper bound of the total mass of neutrinos is £ m v < 0-48 eV. Moreover, 
there exists clear degeneracy between spatial curvature and massive neutrinos in the holographic dark energy model, 
which enlarges the upper bound of £ by more than 2 times. In addition, we demonstrate that, making use of the 
full WMAP data can give better constraints on the holographic dark energy model, compared with the case using the 
WMAP "distance priors". 

PACS numbers: 98.80.-k, 95.36,+x. 
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I. INTRODUCTION 



Observations of type la supernovae (SNIa) [1], cosmic microwave background (CMB) [2] and large scale 
structure (LSS) [3] all indicate that the Universe is undergoing an accelerating expansion. This implies the 
existence of a mysterious component, called dark energy [4], which has negative pressure and takes the 
largest proportion of the total density in the present Universe. In the past fifteen years, lots of efforts [5-13] 
have been made to understand dark energy, yet we still know little about its nature. 

In this paper we focus on the holographic dark energy model, which is a quantum gravity approach to 
the dark energy problem [14]. In this model, the vacuum energy is viewed as dark energy, and is related 
to the event horizon of the Universe when we require that the zero-point energy of the system should not 
exceed the mass of a black hole with the same size [15]. In this way, we have the holographic dark energy 
density [16] 



where c is a dimensionless phenomenological parameter which plays an important role in determining the 
properties of the holographic dark energy, M P \ is the reduced Planck mass, and Rh is the future event horizon 
size of the Universe, defined as 



The holographic dark energy model has been proven to be a competitive and promising dark energy can- 
didate. It can theoretically explain the coincidence problem [16], and is proven to be perturbational stable 
[17]. Moreover, it is favored by the observational data [18]. For more studies on the holographic dark 
energy model, see, e.g., [19-23]. 

It is fairly difficult to calculate the cosmological perturbations of dark energy in the holographic dark 
energy model, because there is a non-local effect making the calculation of perturbations extremely hard 
to be treated. Thus, in the past, for the CMB observations, only the WMAP "distance priors" data were 
used to constrain the holographic dark energy model in order to avoid the inclusion of the perturbations 
of dark energy (see e.g. Refs. [24-28]). However, recently, it has been realized that, in order to make 
progress one should first ignore the non-local effect in the holographic dark energy model and directly 
calculate the perturbations of holographic dark energy as if it is a usual dynamical dark energy. Alone this 
line, in a recent work [29], a global fitting analysis on the holographic dark energy model was performed. 
By using the 557 SNIa data from Union2 sample [30], the baryon acoustic oscillation (BAO) data from 
SDSS DR7 [31], and the full CMB information from WMAP 7-yr observations [32], the author obtained 



p de = 3c 2 M 2 pl R~\ 



(1) 




(2) 
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c = 0.696^-^^^ }^^^ at i ? 2,3cr confidence levels (CL), showing that the future Universe will be 
dominated by phantom dark energy at 3<x CL. However, in Ref. [29] it is not mentioned how the divergence 
problem concerning the w = -1-crossing is treated. 

In this paper, we shall perform a global fit analysis on the holographic dark energy model with spatial 
curvature and massive neutrinos. As argued in Ref. [33], the numerical studies of dynamical dark energy 
should include the current fractional curvature density Q^o as a free parameter to be fitted alongside the 
equation-of-state (EOS) parameters w of dark energy. In addition, the total mass of neutrinos £ m v is a l so 
tightly correlated with w [34]. So, in our work, based on a Markov Chain Monte Carlo (MCMC) global fit 
technique, we will consider spatial curvature and massive neutrinos in the holographic dark energy model, 
and will deeply analyze the influences of these two factors on the fitting results. In a dynamical dark energy 
model, one must be careful about the treatment of perturbations in dark energy when w crosses -1. In this 
work we use the "parameterized post-Friedmann" (PPF) approach [35] implemented in the CAMB code 
following the WMAP team. 

This paper is organized as follows. In Section II, we derive the basic equations for the holographic 
dark energy in a Universe with spatial curvature and massive neutrinos. In Section III, we introduce the 
methodology and the observational data, and then give the data-fitting results. At last, some concluding 
remarks are given in Section IV. In this work, we assume today's scale factor ao = 1, so the redshift z 
satisfies z = a~ l - 1; the subscript "0" always indicates the present value of the corresponding quantity, and 
the unit with c = h = 1 is used. 

II. HOLOGRAPHIC DARK ENERGY MODEL WITH SPATIAL CURVATURE AND MASSIVE NEUTRI- 
NOS 

In this section, we derive the basic equations for the holographic dark energy model with spatial curva- 
ture and massive neutrinos. 

A. Friedmann equations in a non-flat Universe 

In a spatially non-flat Friedmann-Robertson- Walker Universe, the Friedmann equation can be written as 

?>M 2 pl H 2 =p k + Pdm + Pde +Pb+Pv+ Pr, (3) 

where pk = -3M^ Jr is the effective energy density of the curvature component, p^ m , p^, p^, p y and p r 
represent the energy dinsity of dark matter, dark energy, baryon, massive neutrinos and radiation, respec- 
tively. Notice that we adopt the approximate method used in the five-year analysis of WMAP [36]: dividing 
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neutrino component into the relativistic neutrinos and the massive neutrinos, and including the relativis- 
tic neutrinos into radiation component. For convenience, we define the fractional energy densities of the 
various components, 

O - ~ k - Pk O - Pdm O - Pde O - Ph O - Pv O - Pr (A\ 

H z a z p c p c p c p c p c p c 

where p c = 3M 1 pl H 1 is the critical density of the Universe. It is clear that 

a k + a dm + a de + a b + a v + a r = \. (5) 

In addition, the energy conservation equations for the various components take the forms 

p dm + 3Hp dm = 0, (6) 

p de + 3H(p de + p de ) = 0, (7) 

p b + 3Hp b = 0, (8) 

p v + 3Hp v = 0, (9) 

p r + 4H Pr = 0, (10) 

p k + 2Hp k = 0. (11) 



Combining Eqs. (6)— (1 1) together, we can obtain the form of p de , 

2 H 1 

'3lfi Pc ~ Pc ~3 Pr+ 3 



2 # 11 

Pde = --TPJPc ~ Pc~ ~Pr + ~Pk> (12) 



Substituting p de into Eq. (7), we have 

( 2 i + §f + 3H } pde + H(Pk ~ pr) ~ ( 2 i + 3H ^ Pc = °- (13) 

Dividing the above equation by p c , we get a derivative equation of // and Cl de , 

2(Q de - 1)- + Q de + H(3Q de -3 + n k -n r ) = 0. (14) 
ti 
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B. Holographic dark energy in a non-flat Universe 



From the energy density of the holographic dark energy, we have 

L= C -—. (15) 

In a non-flat Universe, the IR cut-off length scale L takes the form 

L = ar{t\ (16) 

and r(t) satisfies 

r® dr r +o ° dt 

; —7. (17) 

Jo Vl - kr 2 Jt a (t) 

By carrying out the integration, we have 

m = ± sin ( V* r *) = -L sin ( V* P -^). (18) 

Equation (16) leads to another equation about r(t), namely, 

L c 

r{t) = -= — . (19) 



Combining Eqs. (18) and (19) yields 

jQ. de aH 

Taking derivative of Eq. (20) with respect to t, one can get 



r- r +co dt C<k 
Vk I — = arcsin — — — . (20) 
Jt a -jQ7 e aH 



^ +H+ */o^p_t (21) 



2^ " H V c 2 a 2 ' 
C. Evolution equations of E(z) and Qde(z) 

Combining Eq. (14) with Eq. (21), we eventually obtain the following two equations governing the dy- 
namical evolution of the holographic dark energy in a Universe with spatial curvature and massive neutrinos, 

1 dE(z) il de 



E(z) dz 1+z 



Q k - Q r - 3 1 / Side ~ 
2Q de + 2 + i~ + ak 



(22) 



dQ.de 2Q. de {\ - Q. de ) [ Q. de _ 1 O k -O r 

~ = U— + ak + 2-20^7) ' (23) 
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where E(z) = H(z)/Ho is the dimensionless Hubble expansion rate, and 

gio(l+ z) 2 Q r0 (l+ z) 4 .... 

^ = -e^-' ^ = (24) 

The initial conditions are £(0) = 1 and Q<fe(0) = 1 - Qyto _ &dm0 - Qbo ~ QvO _ ^rO- Note also that Q v o can 
be expressed as [32] 

- 94^eY' (25) 

where /z is the reduced Hubble constant, and £ m v is the sum of neutrino masses. In addition, the value of 
f2 r o is determined by the WMAP 7-yr observations [32] 

Q r0 = 2.469 x l(T 5 /r 2 (l + 0.2271^//), (26) 

where N e ff = 3.04 is the effective number of neutrino species. Thus, Eqs. (22) and (23) can be solved 
numerically, and will be used in the data analysis procedure. 

In this paper, we shall consider the following four cases: (a) the model of holographic dark energy 
without spatial curvature and massive neutrinos (ClkO = and £ m y = 0), denoted as HDE; (b) the model 
of holographic dark energy with spatial curvature but without massive neutrinos (QkO ^ but X m v = 0), 
denoted as KHDE; (c) the model of holographic dark energy with massive neutrinos but without spatial 
curvature (X m y ± but Q^o = 0), denoted as VHDE; (d) the model of holographic dark energy with spatial 
curvature and massive neutrinos, (Q^o ^ and X m v ^ 0), denoted as KVHDE. 



III. METHODOLOGY, DATA, AND RESULTS 

We have modified the MCMC package "CosmoMC" [37] to perform a global fitting analysis for the 
considered models. As is known, within the dynamical dark energy models, the perturbations of dark energy 
are important [38-40] for the global fitting analysis. In this work we calculate the dark energy perturbations 
by using the formalism of Ma and Bertschinger [41]. Also, we are very careful about the treatment of the 
divergence problem [42] for the dark energy perturbations when w crosses -1. We deal with this issue by 
using the PPF code [35]. This code supports a time-dependent EOS w that is allowed to cross -1 multiple 
times [43] for the dark energy perturbations. Moreover, it has been widely used in the literature to deal with 
the perturbations of dark energy (see e.g. Ref. [44]). 

Our most general parameter space vector is: 

P = (a) b , oj dm , 0, r, c, Q k0 , V ra y , n s ,A s ), (27) 
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where = Clboh 2 and oj^m = Qdm$h 2 , is the ratio (multiplied by 100) of the sound horizon to the angular 
diameter distance at decoupling, r is the optical depth to re-ionization, c is the phenomenological parameter 
of the holographic dark energy model, A s and n s are the amplitude and the spectral index of the primordial 
scalar perturbation power spectrum. For the pivot scale, we set k s o = 0.002Mpc _1 to be consistent with the 
WMAP team [32]. 

In the computation of the CMB anisotropy, we include the WMAP 7-yr temperature and polarization 
power spectra [32] with the routine for computing the likelihood supplied by the WMAP team [45]. For 
the SNIa, we make use of the recently released 580 SNIa data from the "Union2.1" sample [46], where the 
systematic errors of SNIa are included in our analysis. For the LSS information, we use the BAO data from 
the SDSS DR7 [31] and WiggleZ Dark Energy Survey [47]. In addition, we also use the latest Hubble space 
telescope (HST) measurement of the Hubble constant, Ho = 73.8 + 2.4 km s" 1 Mpc" 1 [48]. 

Before constraining the model parameter space, we study the effects of the phenomenological parameter 
c on the CMB Cj T power spectrum by setting different c and fixing the other model parameters. As an 
example, here we just consider the simplest HDE model. The results are shown in Fig. 1. For comparison, 
we also plot the results of the XCDM models with different constant w. One can see that, for the HDE 
models with different c, the main difference appears at low (Z < 20) multipole momentum parts which 
correspond to large scales. This means that the main contribution to CMB power spectra from the HDE 
comes from the evolution of EOS and its domination at the late epoch through late ISW effect. As seen 
in this figure, a smaller c will yield a smaller Cj T at low Z. Besides, since the EOS of HDE satisfies 

w ~ ~5 3c~ I**]' a smaller c will give a smaller w. So for the HDE model, a smaller w will lead to a 

smaller Cj T at low Z. Similarly, for the XCDM model, a smaller w will also lead to a smaller Cj T at low Z. 
Thus, the results of the HDE and the XCDM model are consistent with each other. 

In the following, we will present the data-fitting results. Table I summarizes the fitting results, including 
the best-fit and 1, 2, 3<x values of c and Qko, as well as the 2<x upper bounds of £ for the considered 
models. Moreover, we list the maximal confidence levels for c < 1, in this table. In addition, in Figs. 2-5, 
we plot the ID marginalized distributions of individual parameters, as well as the 2D marginalized 1, 2, 3<x 
CL contours, for these models. Let us discuss them in detail in what follows. 

The fitting results of the HDE model are shown in Fig. 2. We find that the values of c in 1, 2, 3<x regions 
are c = 0.680^ ^^ }^^^^; note that the error bars are slightly smaller than those obtained by using 
the WMAP "distance priors" [28]. Moreover, we find that c < 1 at 4.2cr CL, which means that the future 
Universe will be dominated by phantom dark energy and will end up with a "big rip" (cosmic doomsday) 
at more than Act CL. In Ref. [23] it has been demonstrated that c < 1 may lead to the ruin of the theoretical 
foundation — the effective quantum field theory — of the holographic dark energy scenario. To rescue the 
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FIG. 1 : The effects of the phenomenological parameter c on the CMB Cj T power spectrum. The black dots with error 
bars denote the observed data with their corresponding uncertainties from WMAP 7-yr results. The solid lines denote 
the holographic dark energy models with different c, and the dotted lines denote the XCDM models with different w. 
For the other model parameters, we adopt their best- fit values given by the WMAP 7-yr observations [32]. 



TABLE I: Fitting results of the holographic dark energy models 



Model 


c 


Maximal CL for c < 1 




2 m v (2cr) 


HDE 


fiSO+O.064+0. 135+0.222 
u,uo -0.066-0. 1 19-0. 159 


4.2cr 






KHDE 


r> 7 r>9 +0.104+0.232+0.393 
u ' /u -0.063-0.102-0.176 


2.5cr 


o nnzL + o-oo 9+ o-o i6+ o- 023 

u,uu -0.004-0.010-0.015 




VHDE 


7OQ+0.014+0.111+0.159 
U,/UO -0.099-0.153-0.215 


4.6cr 




2 m v < 0.48 eV 


KVHDE 


O 7QQ+0.037+0.185+0.321 
u ' 1 JJ -0.107-0. 170-0.230 


2.7o- 


n 01 n + o- 01 o + o-o 2 o + o-o 32 

u ' w -0.004-0.014-0.018 


2m y < 1.17 eV 



holographic scenario of dark energy, one may employ the braneworld cosmology and incorporate the extra- 
dimension effects into the holographic theory of dark energy. It has been found [23] that such a mend could 
erase the big-rip singularity and leads to a de Sitter finale for the holographic cosmos. In addition, if there 
is some direct, non-gravitational interaction between holographic dark energy and dark matter, and such an 
interaction satisfies some conditions, the big rip can also be avoided [28, 49]. 

Figure 3 shows the results of the KHDE model. It is seen that, compared with the HDE model, the 
KHDE model slightly favors a larger best-fit value of c. Moreover, the error bars of c are also enlarged. As 
shown in Fig. 2, for the HDE model, we have c < 1 at more than Act CL. However, As shown in Fig. 3, 
after considering spatial curvature, we have c < 1 only in 2.5cr range. 

This figure also shows the degeneracy situation of QkO and c, in the KHDE model. It is clear that QkO 
and c are in positive correlation. This result is well consistent with our pervious work (see Fig. 4 of [49]; 
notice that the convention of Q^o in this paper is different from that in [49]). The best-fit value for Q^o 
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FIG. 2: The ID marginalized distribution of the phenomenological parameter c, for the HDE model. 




FIG. 3: The ID marginalized distributions of individual parameters and 2D marginalized 1, 2, 3<x CL contours, for the 
KHDE model. 

is very close to zero, i.e., QkO = 0.004. The 2<x range of the spatial curvature is -0.006 < QkO < 0.020. 
Actually, the 3cr error bars of QkO are still fairly small, also in order of 10" 2 . 

The constraints on the VHDE model are shown in Fig. 4. One can see that, compared with the HDE 
model, the VHDE model also slightly favors a larger best-fit value of c. However, for the VHDE model, 
the changes on the error bars of c are quite different from the KHDE model: in the KHDE model, both the 
upper and the lower bounds of c are enlarged comparing to the HDE model; while in the VHDE model, 
although the lower bounds of c are enlarged, the upper bounds of c are reduced comparing to the HDE 
model. Moreover, for the VHDE model, we have c < 1 at 4.6<x CL, which is quite similar to the result of 
the simplest HDE model. Therefore, we can conclude that the inclusion of massive neutrinos does not have 
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FIG. 4: The ID marginalized distributions of individual parameters and 2D marginalized 1, 2, 3<x CL contours, for the 
VHDE model. 

significant influence on the phenomenological parameter c. 

From Fig. 4, one can also see that there is no significant correlation between £ m y and c, in the VHDE 
model. In addition, we obtain the upper bound of the total mass of neutrinos in the holographic dark energy 
model, X m v < 0.48 eV at 2<x CL. This is the first result of the neutrino mass in the holographic dark energy 
model. For comparison, we mention here some results of neutrino mass in other dark energy scenarios. For 
a flat ACDM model, i.e., w = -1 and Q*o = 0, Komatsu et al. [32] found that the WMAP+BAO+# limit is 
Yjifiy < 0.58 eV (95% CL). For a constant w model the results given by Komatsu et al. [32] are: £ m v < 0.71 
eV (95% CL) from WMAP+LRG+H , and £ m y < 0.91 eV (95% CL) from WMAP+BAO+SN (where SN 
is the Constitution sample). 

The results for the most sophisticated case, i.e. the KVHDE model, are shown in Fig. 5. We find that, 
for the KVHDE model, we have c < 1 only at the 2.1 cr level. 

This figure also shows the degeneracy situations of various parameters, in the KVHDE models. In the 
KHDE model we find that Q#o and c are in positive correlation, and in the VHDE model we do not observe 
significant correlation between £ and c. So, we believe that the parameters QkO and £ m v should be 
in positive correlation. Indeed, in the KVHDE model, we find that the result is in accordance with the 
expectation. In addition, when simultaneously considering spatial curvature and massive neutrinos in the 
holographic dark energy model, the parameter space of (£^o, 2 m y ) is greatly amplified. For example, the 
upper bound of X is enlarged by more than 2 times comparing to the VHDE model. 

Next, let us discuss the cosmological consequences of introducing spatial curvature and massive neutri- 
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FIG. 5: The ID marginalized distributions of individual parameters and 2D marginalized 1, 2, 3<x CL contours, for the 
KVHDE model. 

nos in the holographic dark energy model. In Fig. 6, we plot the CMB Cj T power spectrum for the HDE, 
KHDE, VHDE and KVHDE models with the corresponding best-fit parameters. To make a comparison, we 
also include the ACDM model with the best-fit parameters given by the same set of data. Two significant 
features can be seen in this figure: first, the Cj T power spectrum for the HDE is well inside the error bars 
of the observational data given by the WMAP 7-yr measurements and matches to the ACDM model very 
well, this implies that the HDE is a competitive model of dark energy; second, the Cf T power spectrum 
for the KHDE, VHDE and KVHDE models almost overlap with the results of the HDE model, showing 
that the inclusion of spatial curvature and massive neutrinos in the holographic dark energy model does not 
significantly change the best-fit results. 

In Fig. 7, we plot the evolution of w(z) along with redshift z, including the best-fit results, as well as 
the 1, 2, 3, Act regions, for the considered models. As seen in the left panels of Fig. 7, for the holographic 
dark energy models without spatial curvature, w will cross -1 at more than Act CL. As mentioned above, 
this means that the future Universe will be dominated by phantom dark energy and will end up with a 
"big rip" singularity at more than Act CL. As seen in the right panels of Fig. 7, after introducing the 
spatial curvature, w will cross -1 only in about 2.5<x range. Therefore, the inclusion of spatial curvature 
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FIG. 6: The CMB Cj T power spectrum for the HDE, KHDE, VHDE and KVHDE models with the corresponding best- 
fit parameter values. The black dots with error bars denote the observed data with their corresponding uncertainties 
from WMAP 7-yr results. 

in the holographic dark energy model may be helpful to alleviate the future cosmic doomsday problem. In 
contrast, the inclusion of massive neutrinos does not have significant influence on the evolution of w(z). 

We are also interested in the differences between the cosmological constraints given by the WMAP 7-yr 
"distance priors" and those given by the full WMAP 7-yr data. As an example, we plot the marginalized 
1, 2, 3<x CL contours in the c-Q^o plane, for the KHDE model, in Fig. 8. The pink contours are plotted 
by using the WMAP "distance priors", and the olive contours are plotted by using the full WMAP data. 
As seen in this figure, the regions of the 1, 2, 3<x contours given by the full WMAP data are significantly 
smaller than the corresponding regions given by the WMAP "distance priors". Therefore, making use of 
the full WMAP data can give better constraints on the holographic dark energy model, compared with the 
case using the WMAP "distance priors". 

IV. CONCLUDING REMARKS 

In this paper we consider the holographic dark energy model with spatial curvature and massive neutri- 
nos. It is well known that both the spatial curvature and neutrino mass are correlated with the dark energy 
EOS, so it is important to study the influences of these factors to the holographic dark energy. In addition, 
it is also rather significant to consider the cosmological perturbations in holographic dark energy and make 
a global fit analysis on the holographic dark energy model. 

We placed constraints on the holographic dark energy in a Universe with spatial curvature and massive 
neutrinos, based on a MCMC global fit technique. The cosmic observational data include the WMAP 
7-yr temperature and polarization data, the SNIa data from Union2.1 sample, the BAO data from SDSS 
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FIG. 7: The evolution of w(z) along with redshift z, including the best-fit results, as well as the l-4<x regions, for the 
considered models. 

DR7 and WiggleZ Dark Energy Survey, and the latest measurements of Ho from HST. In order to treat the 
perturbations in dark energy when w cross -1, we employed the PPF method. So, we do not suffer from the 
divergence problem when w crosses -1. 

We found that, for the simplest HDE model, the phenomenological parameter c < 1 at more than 4<x 
CL, showing that the future Universe will be dominated by phantom dark energy at more than Act CL. 
After taking into account spatial curvature, we have c < 1 only in about 2.5<x range, implying that the 
inclusion of spatial curvature in the holographic dark energy model may be helpful to alleviate the future 
doomsday problem. In contrast, the inclusion of massive neutrinos does not have significant influence on 
the phenomenological parameter c. 

For the KHDE model, we found that the 2cr range of the spatial curvature is -0.006 < QkO < 0.020; 
moreover, the 3<x error bars of Q#o are still fairly small, also in order of 10~ 2 . For the VHDE model, 
we obtained that the 2<x upper bound of the total mass of neutrinos is £ m y < 0.48 eV, which is the first 
result of neutrino mass in the holographic dark energy model. Moreover, when simultaneously considering 
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FIG. 8: The marginalized 1, 2, 3<x CL contours in the c-Q.kO plane, for the KHDE model. The pink contours are plotted 
by using the WMAP7 distance priors, and the olive contours are plotted by using the full WMAP7 data. 

spatial curvature and massive neutrinos, the upper bound of £ will be enlarged by more than 2 times. 
Furthermore, we also demonstrated that, making use of the full WMAP 7-yr data can give better constraints 
on the holographic dark energy model, compared with the case using the WMAP "distance priors". 

It should be mentioned that, there are still some factors not covered in our paper, e.g., the interaction 
between dark sectors. Evidently, when taking into account the interaction, the computation of the dark 
energy perturbations will become much more complicated. These issues deserve further investigations in 
the future work. 
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